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Abstract 

Various ensembles of random matrices with independent entries are analyzed by 
the replica formalism in the large- limit. A result on the Laplacian random matrix 
with Wigner-rescaling is generalized to arbitrary probability distribution. 



1 Introduction 

Since the introduction of bosonic replicas in the study of disordered systems and ran- 
dom matrices several decades ago [1] , most analytic investigations followed a path with 
three basic steps : 

1) After defining a replica partition function, it is easy to average it over the ensemble 
of random matrices to obtain a replica partition function where the interaction of the 
replica fields depends on the probability distribution of the ensemble. 

2) A set of auxiliary fields are introduced by Hubbard- Stratonovich identities to decou- 
ple the interacting pairs, transforming the replica partition function in a form suitable 
to the thermodynamic limit. In recent years a Gaussian functional identity conve- 
niently replaced the set of auxiliary fields. 

3) The saddle-point method leads to a non-linear integral equation of very difi'erent 
complexity : it is trivially solved in the case of Wigner matrices whereas it is extremely 
difficult in the case of the Laplacian of a random graph. 

This paper follows the well established steps, respectively in sections 2, 3, 4, with 
greater generality. Instead of choosing a definite probability distribution, we consider 
the off-diagonal matrix elements to be independent identically distributed (i.i.d.) ran- 
dom variables, where their probability distribution is coded by the cumulants. The 
diagonal matrix elements will be i.i.d. random variables statistically independent from 
the off-diagonal matrix elements or, as in the case of the Laplacian matrix, statistically 
dependent on them. In either case (statistically dependent or not) different specifi- 
cations of the cumulants correspond to a full matrix or a sparse one. This general 
approach does not lead to complications and all the three above mentioned steps are 
performed. 

In sect. 3 we discuss some puzzling features of the useful functional representation that 
decouples pairs of replica fields and its equivalence with the older method of auxiliary 
variables. The reader who is satisfied with formal inversion of operators without closer 
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inspection, may safely skip reading this section. 

Sect. 4 contains the saddle-point method , the non-linear integral equation and the new 
result of the paper : the limiting resolvent of the Laplacian of a random matrix with 
i.i.d. entries, rescaled a-la-Wigner, may be computed by a naive use of Pastur 
law of addition of random matrices, despite the statistically dependence of the diago- 
nal entries. Our new result generalizes a previous one [2], obtained by diagrammatic 
methods and the assumption that the off-diagonal matrix entries are i.i.d. central 
Gaussian variables. 

Wc think that the general and straightforward approach of this paper has some 
merits : indeed, as a by product, we obtain, at the beginning of Sect.4 the well known 
Pastur law of addition of random matrices, in the replica formalism. 

Sect.4 ends with the non- linear integral equation pertinent to the random graph. 
We included it to emphasize the generality of the replica formalism and the difference 
between Wigner-rescaled matrix ensembles (and their Laplacians) which leads to ana- 
lytic solutions whereas ensembles where the matrix entries are of order 1/N are much 
more difficult. 

It seems likely that everything in the paper may equally be expressed in terms of 
the supersymmetric method. 



2 Exact representations 

Notation. 

Let H he a, N X N real symmetric matrix. Let us suppose the diagonal entries Hrr are 
TV i.i.d. random variables and the probability density of any of them is Pd{Hrr)- Let 
us suppose the off-diagonal entries Hr^s, where r < s, arc {N^ - N)/2 i.i.d. random 
variables and the probability density of any of them is Poff{Hr,s). Let us write the 
expectations 

< d™ >= j {Hrrr Pd{Hrr) dHrr , < v"" >= j (i/..)'" PoffiHrs) dHrs 

Let G{z) be the matrix resolvent 

Gij{z) = {{zlN-H)-'').. , lmz>0 (2.1) 

Perturbation expansion in powers of H shows that the expectations of all diagonal 
elements < Gii{z) > are equal and of the ofF-diagonal elements < Gij{z) > are also 
equal. 

<G,,(^)>=| V-^^; ' j,<trGiz)>=g{z) 

^ ' \ giiz) if J 

The representation of the inverse matrix and usual definitions for the partition function 
Z{z) and the averaged spectral density p(A) are 
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^ij — 



N/2 



y/det{zlN - H) 

g{z) = -j^-^<\ogZ{z)> , p(A) = -ilmff(A) (2.3) 



To get rid of the denominator of the resolvent, let us introduce n copies 
a = 1,2, .., n 



Z„{z) = Z"{z) 



inn- 

-oo \ 1 



, s 2 8..^ Zr,{z) - 1 

N az n-.o n 



(2.4) 



Exact representations. 

The above introduction of bosonic replicas allows an easy averaging over the matrix 

ensemble and leads us to an exact representation for any N for the averaged < Zn{z) > 
and resolvent, see eq. (2.7) below, in terms of cumulants of the probability distribution 
of the matrix entries. 



Let us assume we may first average on the entries Hij , then integrate over 



finally take the limit n — » 0. Let us split the diagonal from the off-diagonal entries 



n 



< e 



•E„ 



.(o) ,(o) 



> 



(2.5) 



This has a compact representation in terms of cumulants of the moments of the 
matrix entries. Let a; be a random variable, we may define the cumulants Cj by equating 
the powers of a of the exponential generating functions 



log Ke"" >= log 



k 

EC^ k 



that is 



Cl =< X > , C2 =< > — < X >^ 



C3 =< x^ > -3 < x'^ X X > +2 < X >^ 



C4 =< > -4 < >< X > -3 < x"^ >^ +12 < x'^ X X >^ -6 < X 



To simplify the notation, we call c^^d, Ck.v the k-cumulants of the diagonal and 
off-diagonal entries of the random matrix H. 



(2.6) 
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One obtains the exact representation ^ 



l<r-<s<JV 



(2.7) 



The random Laplacian matrix 

An analogous integral representation holds for a random Laplacian matrix. 

Let us suppose that the real symmetric random matrix Hij is a Laplacian matrix: we 

assume that the diagonal elements are Hu =^ — X^^-^j Hij. Again we assume that the 

A'^(A'^ — l)/2 random variables Hij with i < j are i.i.d. 

Since 

<t>iHij<t>j = - ^(<^i)' ( ^ Hij + YHij\+2Y^ <i>i4>jHij = 

= (2.8) 



the partition function , the averaged resolvent and their replicated versions are 



Z„iz) 



\a=li=l / 

,(.) = - 2 ^H^<M£)^> (2.9) 
N dz n^o n 

Again one may perform the average over the matrix ensemble to obtain an exact 
expression in terms of the cumulants Cj of the moments of the off-diagonal entries 



E E' 

l<r<s<N j=l ' 



E' 



(2.10) 



^It is easy to check that this integral representation is an absolute convergent in- 

, J ~ 



tegral for every integer n. Indeed 



< e 



< 1 and similar relation holds for the factor containing Hr,. 
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Special cases. 

It may be useful to recover from the general representations in eqs.(2.7) and (2.10) some 
cases corresponding to special probability distribution for random matrices, which have 
been studied by use of bosonic replicas. 

1. The Gaussian ensemble. In the Gaussian orthogonal ensemble (GOE) the joint 
probability density of the real symmetric random matrix is P{H) 

P{H) = cons. X e~2^*^''^ 

Then all cumulants are zero except the second one, C2,d = and C2,v = cr^/2. 
The exponent F in eq.(2.7) is 



(2.11) 



The F exponent of the replicated partition function of the Laplacian matrix with 
Gaussian entries (2.10) is 



=-- y 

16 ^ 



(2.12) 



2. The sparse random matrix. In the case of the sparse random matrix considered 

by Rodgers and Bray [3] , the independent entries Jij, i < j, have a oven probability 



distribution so that < v 



>= 0, < v^" >= p/N for fe > 1. The same even 



probability distribution and its asymmetric generalization have been studied in [4]. 
The equation to evaluate the cumulants is 



log [< ( 



log 



l + |:(-l + cosa)] =^|(-ia)^- 



Then the exponent F in eq.(2.7) is 



+ 



^10, 

r 

l<r<s<Af 



P_ 

N 



- 1 + cos 



iE( 



l + |p -l + cos5^,/.("Vi"^ 



(2.13) 



If only the limit A?^ — »■ oo is investigated, the first term in the above equation may be 
dropped and the second one may be approximated 



F-l^log 



l + ^(_l + cos^. 



P 

2N 



(2.14) 



3. The random graph. The case of the random graph was often analyzed as a 
Bernoulli probability distribution of the off-diagonal entries of the random matrix 
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Hij, i < j, where Hij is — 1/g with probabihty q/N and with probabihty 1 — q/N , 
while the diagonal entries vanish , Hjj — 0. 

Then < v*" >= (— 1/g)* for > 1 and the equation to evaluate the cumulants cj is 



log < e"^ >= log 



Then the exponent F in eq.(2.7) is 



(2.15) 



(2.16) 



The F exponent of the rephcated partition function of the Laplacian matrix (2.10) is 



F= J2 log 



(2.17) 



l<r<s<N 

The spectral density corresponding to eq.(2.17) is studied in [5], [6], [7]. 



4. The diagonal random matrix. For future reference we also recover the partition 
function and resolvent for the simple case of the diagonal random matrix. Of course 
in this trivial case, replicas are totally unnecessary. 



Prom eq.(2.6) and (2.7) 



. a d ^ 

< e >= 



E^r<<''>= 



fc=0 



< Z„{z) > 



i,oc r=l k=0 

/n oc 
a=l k=0 

.k=0 oc=l 

oo 

p')"'TM<''*>i(i+>)-(i+'-' 



^in7r/4 

Then for n -> 

< Zn{z) >- 



(2.18) 



> z 



k=l 



z-\ 



(2.19) 
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3 Hilbert- Schmidt with repUcas 



The next step in the analysis of the replicated partition function , eqs.(2.7) and (2.10), 
is to transform the part of the integrand of < Zn{z) > that involves sums on pairs of 

sites into a form involving only sums over single sites. Having done that, the partition 
function may be written as an integral where the rmmber of sites A'^ appears just as a 
parameter. 

This was done by Rodgers and Bray [3] by introducing a set of auxiliary variables 
<la,Qaf3,qa0-y, ■■ with the corresponding Hubbard-Stratonovich identities. The same 

authors used a functional representation in the case of a Laplacian matrix, [8]. 
It seems very convenient to use a functional Hubbard-Stratonovich transform 



suggested by Fyodorov and Mirlin [9] and discussed in more detail in the appendix of 
[10]. 

An analogous functional integral representation is discussed by Parisi [11] for euclidean 
random matrices 

e-^E.,3^(-'— ^) = (consi.) j rf^e^/''-<^^-W'^-'(-^)-W+E,-(-^) (3.2) 
As it is indicated by Parisi, eq.(3.2) follows from the familiar identity 

(const) j da j'dxdy(y(x)V~^(x-y)rj{y) + j'dxJ(x)rj{x) _ ^-/S j dx dy J (x)V (x-y)J (y) 

(3.3) 

after the choice of the source J{x) = X^^j^ '^(a; — Xj). 
The problem 

In the Gaussian case, the sum over pairs of sites, occurring in the exponent F, eqs.(2.7) 
and (2.10), is a homogeneous polynomial of degree 4. In the case of a general prob- 
ability distribution, cumulants of arbitrary high order are not zero. Yet it will be 
apparent in the next section that also in this case the large- iV saddle point analysis of 
the resolvent leads to the Gaussian case in leading order : cumulants Cj,^ with j > S 
are associated to non-leading terms, if the random matrix elements are rescaled a-la- 
Wigner. 

Then these integral representations are rather puzzling because they involve the in- 
verse of the integral operator K{4>r, <t>a) which is a polynomial (in the most important 
case) and therefore the eigenspace corresponding to the null eigenvalue has infinite 
dimension. 

Still the method is very useful and leads to correct results. It was correctly indicated 
by Fyodorov how to evaluate the inverse operator in a subspace. This section is a 
little elaboration on Fyodorov appendix [10]. 

Our conclusion is that, provided a regularization of the problem is allowed, the func- 
tional Hubbard-Stratonovich transform is correct and identical to the set of usual 
Hubbard-Stratonovich identities. 
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The symmetric matrix 

Let us consider all replicas in a box , -L < 0^"^ <L,r = 1,..,N, a = 1, .., n. 
Let $r be the n-component vector with components </>r™\ a = 1, .., n. 
Consider the integral operator kernel K{(j)r,4>s) with finite rank D (to make the dis- 
cussion simpler, we shall ignore the imaginary unit factors in this section) 

D / n 

K[i, = E I E = E I (^^- = 

i=\ \a=l / j=l 

3=1 ai,..,oij 

D / N \ 2 

E^(0.,0.)=E| E yL^'r'---^'^''] (3.4) 

r,s j = l f^lj'-iOij \r=l / 

It maps an arbitrary integrable function f{4>a) to the function g{4>r) € So which is the 
linear finite Z)-dimensional vector space spanned by the set s of linearly independent 
vectors 

/L i-L n 

.. K{$r,$s)f(.$s)d$s=g(.$r) , =JJ#<"' (3.5) 

-L J-L ^ 

The vector bases s is made of D non-orthogonal vectors bu(4>r) ■ The integral operator 
K{<pr,(j)s) has D orthonormal eigenvectors Cu, corresponding to non-zero eigenvalues 

L 

K{4>r, ct>s) e.v{4>s) d(j>s = \vev{4>r) , 

■L 

D 



1^=1 fl=l ^ 

■'-^ i/=l 

(3.6) 

The vector bases {6i/(0r)} , and the orthonormal bases {eu{4'r)\ are related by a 
linear invertible matrix A 

bi, = '^A^^eij, , ei, = E(^~^)i, ''m > Avij. = / 6^(0r)e^(^T-)#r- (3.7) 

"'^ J-L 

MM 

As indicated by Fyodorov, the functional Hubbard- Stratonovich transform (3.1) is 
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a compact representation of the multiple integral 

D r_«, . -,2 D 



After a trivial relabelling of the coefficients in the form (3.4), one uses the usual 
set of Hubbard-Stratonovich identities 



1 J —oo V27r 



Finally one may notice that the set of integration variables Qi, and z^y arc related 
by an orthogonal matrix A which shows the direct equivalence of the two multiple 
integrals (3.9) and (3.10) 



M=l 

D 



a -V^ z A -J^A 



^ArAcT^^ra , X/ ^ X/ ' = ^'^^ (3.11) 

The orthogonality of the matrix A (last line above) follows from the equation 

Kiar = ^ c{v)A^rA„a (3.12) 

which is obtained from 

K{$r,$s) = ^ A,.e^(0r)e^(0s) = ^c(z/)6,.(0r)6^(0s) (3.13) 

after multiplying and integrating J d(f)rd4>s ea{4>r) eT{(t>s) 
Finally 



9{4>r) = s/X,Zueu{4>r) = a/aT [(-4. ^)i,<Tg<T] [(^ ^)i/M^M(<?r)] = 

= ^ -\/Ar-4o-''9o-(^~^)i'M&M(0'-) = y^.Vcjlqi^bf.i^r) (3.14) 



^In the functional integrals, T^g refers to the product, properly normalized, of the coordi- 
nates of the function g{4') with respect to an orthogonal bases, then "Dg = ^ Xv/2-irdzi, 
and eq.(3.1) should be written 



(3.10) 
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The random Laplacian matrix. 

The replicated partition function for tlie symmetric random Laplacian matrix is (2.10). 
Again ignoring imaginary unit factors, the finite rank kernel of the integral operator 

is K{(j)r, <f)s ) 

D 

Let us assume the box rcgularizatioii, as in the discussion of the symmetric random 
matrix above. Just as before, the kernel K{(t>r, 4>s) in eq. (3.15) is a symmetric bounded 
integral operator that maps the integrable function f{(j}a) to the function g{4>r) £ Sd, 
spanned by the non-orthogonal bases b„{(f)r), see eq. (3.5). 

It is then possible to represent the kernel K{(f>r, 4>s) in eq. (3.15) with its orthonormal 
bases and obtain Fyodorov functional Hubbard-Stratonovich representation or repre- 
sent it on the non-orthogonal bases hu{4>r) , perform the sum over couples of sites and 
use the Hubbard-Stratonovich identities. 

Let us consider the first term in the sum in eq.(3.15). By the trivial identity 

(pr ■ (pr + >f>s ■ (l>s = (l -h ^T- • (pr) [l + (ps ■ <Ps) — {(pr ■ (pr) {<Ps • (ps) — 1 

The first term is rewritten as sum of factorized and symmetrical terms 

Y {<fr- (ps) ■ {(pr - 4>s) = 

= y { (l + ■ 0r) (l + ■ $s) - {$r ■ ^r) {^s • (^s ) " 1 " 2 • 0^) } (3.16) 

All further terms in the sum are integer powers of the first term, then a finite sum of 
factorized and symmetrical terms is obtained. 



In particular, in the Gaussian case, where only the cumulant C2 7^ 0, we have 



'^K((pr,(ps) = -^'^[{(pr - (ps) ■ {(pr - (ps)] = 
r,s r,s 

^ {$r ■ $r + $s ■ $s) +4 {$r • <As) — 2 {^r ' 0s) {^r • 0r + 0s • 0s) 



C2 



-2 



^(1 + $r ■ $r) 
y^^$r- (f 



- 2 



^(0r • 0r-)' +N'-2 

r 

^^(0r • 0r)(l + 0r • 0r) 



^ (l 0r • $r) 



a/3 



+^E 



■^E 



E' 



This form is suitable to linearize the squared sums of single sites by the introduction of 
9 auxiliary variables through Hubbard- Stratonovic identities. The method also works 
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for arbitrary power j > 2 of the basic form in eq.(3.16), but it is cumbersome and by 
far less efficient than the functional integral representation. 



4 Large N 

As a preparation for the new result of this paper, about the limiting resolvent of 
Wigner-rescaled random Laplaciams, we first consider an ensemble of Wigner-rescaled 
matrices which leads to a derivation of Pastur addition law of random matrices in the 
replica formalism. 



Pastur addition law. 

The first problem here considered is the large N behaviour for the symmetric matrix 

where the probability distribution of the off-diagonal entries is central Gaussian with 
Wigner scaling, while the probability distribution of the diagonal entries is arbitrary 
and not rescaled. Then Cj^v = iff j 7^ 2 and C2,v — » C2/N. 

After separating the sums over single sites from sums over couples of sites, the repli- 
cated partition function (2.7) becomes 



I 



(4.1) 



By using the functional representation (3.8) we have 



I[9] 



def 



d(f) e 



def C2 

~ 2N 



(0r • $s) 



(4.2) 



The large N behaviour of the functional integral is dominated by the saddle point 
function gs.p. 



Qs.p. {ij) = - 



N 



I[9s.v] 
iN 



l[gs.v. 



3=as-p- 



The above equation shows that gs.p.{ip) is a quadratic function of il). We consider the 
ansatz symmetric in the replica fields 



(4.3) 



9s.p.{<t>) = -f3{z,Cj,d)4>- 



(4.4) 
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then the function P{z,Cj^d) is solution of the nonlinear equation 



feii)=_£aJ_|.7(;3) where 

m = I d0e^<^-«^'-^^+S.=i (4.5) 

This leads to the evaluation of the limiting resolvent because 

Jim = -2 lim = l/3(.,c,,.) (4.6) 

and, by comparing J(/3) in cq.(4.5) with the partition function < Zn{z) > of an 
ensemble of diagonal random matrices in eqs.(2.18), (2.19), we see that ^ 



7(/3) = [< Z„{z - 13) >f^^ , lim g{z) = lim Ga {z - C2g{z)) (4.8) 

JV— >oo N—i-oo 

where Gd{z) is the resolvent of an ensemble of diagonal random matrices with proba- 
bility distribution coded by the cumulants Cj^d- 

The functional equation for the limiting resolvents (4.8) is the law of addition of ran- 
dom matrices proved by L. Pastur [12] 

,(.) = /.A_^^ , MX)'=' -'-ImGdiX) (4.9) 

We remark that the result (4.9) holds unchanged if one considers the general case 
of non-Gaussian probability distribution for the off-diagonal matrix entries, provided 
Wigner scaling is still considered and the first cumulant vanishes, ci,^ = 0. Then 
Cj,v/N^^^ for j > 2. Indeed, for a general probability distribution, the previous 
analysis is modified by the replacement 

oo 

r,s r,s J=2 r,s 

The series of new terms, inserted in eq.(4.3) only add contributions to the limiting 
resolvent, which are non-leading order in the large- A'^ limit. 

The Wigner-rescaled random Laplacian matrices. 

The large- Af spectral distribution of a real symmetric Laplacian random matrix, where 
the matrix entries are Wigner rescaled, may be determined by the same saddle point 
method described above. 



■^If all diagonal cumulants vanish, cj^d = 0, the Gaussian integrals at the right side of 
eq.(4.5) are easily evaluated and one obtains the quadratic equation 0^{z) — z/3{z) + C2 = 0. 
Together with eq.(4.6), one then recovers the semicircle law 

lim g{z) = -p{z, Cj,d) = :^(^- V^^ - 4c2) (4.7) 

JV— ►oo C2 2C2 \ J 
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Let us first consider the simple case of Gaussian probability distribution. Then 
eq.(2.10) becomes 

< Zn{z) >= J e^'^'-^^'-^'-^-^^'-Jif'-^^yi^'-^')]" (4.11) 
Instead of eq.(4.2), we now have 

<Zniz)> = VditCy"l?(7e-5(9.C9)+JVlog/[g] 

K{$rjs) = ^[{$r-$s)-{$r-$s)]' (4.12) 

The large iV behaviour of the functional integral is dominated by the saddle point 
function g^.p. 



I\9s 



SgW 



g=ge.p. 



ii^gi^'?-i?-iSs.p.(i?) 



gs.Ai^) = I [(0- VO • ei^*^-^'---"'^^) (4.13) 

The above equation shows that gs.p.ii') is a quartic function of ip. We consider the 
ansatz, symmetric in the replica fields, 

gs.p.iif) = a{z) + ^f3{z)^f- ■0 + 7(2) (1? • 1?)^ 
then 7(2) = —102/8 and the function f3{z) is the solution of the nonlinear equation 



/(/3) = / # e^^"-")^-^-^!^-^ (4.14) 



This leads to the evaluation of the limiting resolvent because eq.(4.6) still holds and 
lim g{z) = lim Gd,g {z - C2g{z)) (4-15) 

where Gd,g{z) is the resolvent of an ensemble of diagonal random matrices where the 
entries are i.i.d. centered Gaussian variables with variance C2. 

The result (4.15) was derived in [2] by diagrammatic technique and Wick's theorem, 
to account for the Gaussian matrix entries with a constraint on each row of the matrix. 



The present derivation shows that the result (4.15) for the limiting resolvent holds 

also in the general case where the off-diagonal entries of the real symmetric random 
Laplacian matrix are i.i.d. random variables with arbitrary probability distribution, 
coded by the cumulants Cj,„, provided they are Wigner rescaled and that the first 
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cumulant vanishes, ci,v = 0. 

Indeed, for a general probability distribution, the previous analysis is modified by the 
replacement 

K {$r, = X] ^ [(^'' ~ ^'') ' (^'^ ~ ^^)] ^ ~^ 

r,s r,s 
oo 

^ - E ^ E [I (^^^ - ■ (^^^ - (4.16) 

The series of new terms, inserted in eq.(4.13) only add contributions to the limiting 
resolvent, which are non-leading order in the large- TV limit. 

Finally we remark that our result (4.8) refers to an ensemble where the diag- 
onal elements are statistically independent from the off-diagonal entries. Then we 
merely derived, by a different technique a specific case of Pastur's law of addition 

of random matrices. The similar result about Laplacian matrices (4.15) is different 
because, at any finite A'', the diagonal matrices are not statistically independent from 
the off-diagonal entries and Pastur theorem does not apply. Only at iV = oo, because 
of the central limit theorem, the diagonal entries become centered Gaussian random 
variables. Our saddle point analysis indicates that the rate of convergence for the 
diagonal entries to become statistically independent from the entries on the same row 
is fast enough to justify a naive use of Pastur's law also in this case. 

Let L be a real symmetric Laplacian random matrix where the off-diagonal matrix 
entries are i.i.d. random variables. The graphic-combinatorial method described in 
[13] may be used to evaluate < TrL*' >, for any N and generic probability distri- 
bution. We evaluated the Taylor expansion of the resolvent g{z) in 1/z up to order 
(l/z)^. Next it is trivial to rescale the results according to Wigner rescaling of the 
matrix entries, evaluate the large- A'^ limit and find results in agreement with eq.(4.15). 
This adds independent support to the saddle point analysis of this paper. 

The random graph 

In models of sparse matrices or in the random graph model, the probability distribution 
for the entries of the random matrix P{Hij) is such that the cumulants Cj are of order 
1/N , in the large- A?" limit, for every j. The general saddle point method of this section 
is still valid, but the non-linear integral equation obtained for gs.pXi') is more difficult. 
In leading order the diagonal cumulants Cj,d may be neglected and instead of eq.(4.3) 
we have 



I[gs.p]= y"d<^e5"*"-«'-'S''fW (4.17) 

If one looks for solutions with spherical symmetry in the replica space, that is 
gs.p.{i>) = S's.p.dV'Di the angular integrations are easily done and only the terms in the 
series (4.17) which are proportional to Cj,„ with j even integer contribute. We find 



j=i 
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zp -igs.p.ip) 



p oo / poo 

9s.Am = j p"-'ei^''"-">-''-^^^F{p\i;\)dp/ J p"-'ei 

OO 

fc=l ^ ' 

n l-3---(2fc-3)(2fc-l) 
* n(n + 2) ••• (n + 2/c - 4)(n + 2/c - 2) 

For n~0 T p""^ e*-^"-- ^ ~ i a.. / 1 \ 3 ■ 5 • • • (2fc - 3)(2fc - 1) 

° " ^ n ' UJ 2.4...(2fc-4)(2fc-2) 



and we obtain a simpler non-linear equation 



4 



/•oo 

oc 

/ (^1^1) = E (fc - l)?fc!2^^-i (-1)"^^"^?" (4.19) 
fc=i 

A variety of approximate solutions may be found in [3] and [4]. An analogous 
non-linear integral equation is obtained for the random Laplacian. Its approximate 
solutions may be found in [8] , [6], [7], [5]. A numerical solution is studied in [14]. 

5 Conclusion 

The limiting spectral distribution for ensembles of real symmetric random Laplacian 
matrices with independent identically distributed Wigner-rescaled matrix entries has 

been derived in this paper. This is a generalization of previous results obtained with 
different techniques in the case of Gaussian matrix entries. 

The replica method, combined with a Hubbard- Stratonovich functional identity ( 
reviewed and explained in Sect. 3) turns out to be a very powerful and general method 

in the study of random matrices. In the derivation of our new result, the Pastur law 
of addition of random matrices is recovered by the replica method. 

Our derivation, up to the large- iV saddle point analysis, applies to Wigner rescaled 
matrix ensembles as well as to ensembles of sparse matrices , where usually < >~ 

0(1/N). The relevance of the different scalings emerges at that point: in the case of 
ensembles of sparse matrices, the non-linear integral equation is very difficult to solve 



^As a simple check of our evaluation (4.19), one may insert c^k = p/N + 0{1/N'^) and 
obtain / (pl'i/'l^ = —{p/N)\'ip\Ji{p\'tp\) which reproduces the non-linear equation first derived 
by Rodgers and Bray [3]. 
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analytically. 
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